On Breaking Time Reversal in a Simple, Smooth, Chaotic System 
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Within random matrix theory, the statistics of the eigensolutions depend fundamentally on the 
presence (or absence) of time reversal symmetry. Accepting the Bohigas-Giannoni-Schmit con- 
jecture, this statement extends to quantum systems with chaotic classical analogs. For practical 
reasons, much of the supporting numerical studies of symmetry breaking have been done with bil- 
liards or maps, and little with simple, smooth systems. There are two main difficulties in attempting 
to break time reversal invariance in a continuous time system with a smooth potential. The first 
is avoiding false time reversal breaking. The second is locating a parameter regime in which the 
symmetry breaking is strong enough to transition the fluctuation properties fully toward the broken 
symmetry case, and yet remain weak enough so as not to regularize the dynamics sufficiently that 
the system is no longer chaotic. We give an example of a system of two coupled quartic oscilla- 
tors whose energy level statistics closely match those of the Gaussian unitary ensemble, and which 
possesses only a minor proportion of regular motion in its phase space. 



Since its introduction into nuclear physics by Wigneni, 
random matrix theory (RMT) has grown to encompass 
a broad variety of applications, and can be viewed as 
a significant portion of the foundation of the statistical 
mechanics of finite systems. One of the central tenets 
of RMT emphasized by Wigner was that the ensembles 
carry no information other than that required by the 
symmetries of the system. A very important symmetry is 
time reversal which determines whether the ensemble is 
composed of real symmetric matrices, Gaussian orthogo- 
nal ensemble (GOE - good symmetry), or complex hermi- 
tian matrices, Gaussian unitary ensemble (GUE - broken 
symmetry). This distinction leads to quite different pre- 
dictions for the statistical properties of both the energy 
levels and the eigenfunctions. In fact, Wigner proposed it 
as a test of time reversal invariance in the strong interac- 
tion using slow neutron resonance dataS. His suggestion 
was not fully realized until more than twenty years later—. 

Time reversal is well known to be an antiunitary sym- 
metry in quantum mechanics, and Robnik and Berry gen- 
eralized the criterium for expecting the statistical prop- 
erties of the GOE to include any antiunitary symme- 
try. An elementary example of an antiunitary symmetry 
would be the product of some reflection symmetry and 
time reversal. Some of the early investigations of non- 
invariant systems were of Aharanov-Bohm chaotic bil- 
liards^, symmetry breaking quantum maps«V), and com- 
binations of magnetic and scalar forcesi 

There is no known mechanical- type system of a particle 
moving under the influence of a simple, closed, smooth 
potential whose dynamics is rigorously proven to be fully 
chaotic, independently of the question of anti-unitary 
symmetry; we exclude diffusive dynamics associated with 
random potentials. For some period of time, the x 2 y 2 
quartic potential was a prime candidate, but eventually 
stable trajectories were found£. Nevertheless, there ex- 
ists a family of quartic potentials with x 2 y 2 as a lim- 
iting case whose fluctuations closely approximate GOE 
statistics, and whose classical dynamics contain negli- 



gible phase space zones of stable trajectories. There- 
fore, to construct a close approximation of a Hamilto- 
nian without an anti-unitary symmetry and GUE statis- 
tics, we begin by considering the symmetry preseverving 
two-degree-of-freedom coupled quartic oscillators whose 
Schrodinger equation is given by 
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-#(r) + V(r)^(f) = E^{r) (1) 



The potential can be expressed as 



V{r) = o(A) 
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where a(A) is a convenient constant, b ^ 1 lowers the 
symmetry, and A gives the strength of the coupling; 
(b = 1, A = —1) is equivalent to V(r) — x 2 y 2 by a 7r/4 ro- 
tation of the [x, y) coordinates. The corresponding clas- 
sical Hamiltonian is 
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H (r,p) = ^- + V(r) 
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This system is symmetric under time reversal and re- 
flections with respect to x and y. For strong couplings, 
A < —0.6, the statistics have been shown to agree ex- 
tremely well with the GOE predictions^. As a first at- 
tempt to break time reversal symmetry, we add the fol- 
lowing term to the quantum Hamiltonian (H — Hq+Hi), 
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where r is the radial polar variable. Hi breaks time re- 
versal invariance in the Hamiltonian without altering the 
original reflection symmetries, and thus does not admit 
an antiunitary symmetry from any combination of reflec- 
tion and time reversal. This term was chosen to maintain 
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a scaling property of the eigensolutions of the quartic os- 
cillators due to the homogeneity of the potential 9 . How- 
ever, Hi gives us an excellent example of false symmetry 
breaking. In fact, it turns out that Hi is given by the 
cross terms arising from a vector potential A(r) = er 2 f 
(which can also be expressed as the gradient of a scalar 
function); i.e. 

F 1 = ^(V.A(f) + A(r).V) (5) 
Im 

This means that, up to the addition to the potential of 
a term e 2 r 4 , Hi can be understood as deriving from a 
magnetic field B = VxA = 0, which obviously will not 
change any physical quantity. And indeed, it is possible 
to make a gauge transformation and rewrite the wave 
function using the Dirac substitution 



ty(r) = exp 




A(r") • dr* *'(r) 
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in order to cancel Hi plus the aforementioned e 2 r term. 
Furthermore, the transformation is single valued since 
B = implies 



j>A{?) -df 



= 



(7) 



The expectation, following the considerations o^, is 
therefore to find GOE statistics and not those of the 
GUE. See the upper panel of Fig. Q, which compares 
the number variance of the system (i.e. the variance of 
the number of levels found in an energy interval of width 
s scaled locally to mean unit level density) with the pre- 
dictions of the GOE and GUE. The GOE statistics are 
closely matched to a mean level spacing and a bit be- 
yond. The parameter e was chosen slightly greater than 
unity because this forced the eigenlevels through a couple 
of avoided crossings, which would have been sufficient to 
push the spectral fluctuations toward GUE if the sym- 
metry was not being falsely broken. 

From a classical perspective, the Hamiltonian, not in- 
corporating the jauge transformation, is 
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FIG. 1: The number variance for the quartic oscillators 
(solid lines) compared to the GOE and GUE (dashed lines). 
The upper panel shows the statistics for the hidden symme- 
try case, i.e. with the Hamiltonian term in Eq. (4)) with 
(e = 0.5, A = —0.65). The lower panel shows the results using 
the Hamiltonian term of Eq. (11) and (e = 1.0, A = —0.80). 
Spectra containing a total of 600 levels were used to generate 
the statistics. The lowest 50 eigenvalues were dropped. 



H(r,p) 



p 2 - 2er 2 p r 
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+ V(r, 6) (8) 



which itself appears to violate time reversal symmetry, 
as well as in Hamilton's equations of motion: 
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However even without considering a canonical transfor- 
mation, just by converting to the Lagrangian description 



of the dynamics using the left hand side equations, it 
turns out that 



L = —Imr 
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)-V(r,9) + —r 4 + er 2 r (10) 
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The final term, which appears to break the symmetry, 
cannot enter the equations of motion. They are invariant 
under addition of any total time derivative. It is not 
always obvious, a priori, whether a symmetry breaking 
term leads to false symmetry breaking or not. 

If we multiply the symmetry breaking term by any 
function of 9, it can no longer be a total time derivative, 
nor can the vector potential be expressed as the gradient 
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(X = -0.8 , e = 1.0) 
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FIG. 2: The x — Poincare surface of section for the 
Hamiltonian with the term from Eq. (11), and parameters 
(A = — 0.8, e = 1) as in the bottom panel of Fig. (1). The 
coordinates y and p y are in scaled unit. 



of a scalar function. Consider 

- , ihe 2 . / d d n\ . . 



The vector potential becomes A(r) = er 2 cos 2 Of, and 

therefore B = 2exy/r ^ 0. The integrals J r A(r f ) ■ dr* 
and f A(r*) • dr" are path dependent. The Dirac substi- 
tution is not useful, and false symmetry breaking is not 
an issue. 

Driving the eigenvalues at numerically attainable en- 
ergies through a couple of avoided crossings forces e to 
be chosen in the neighborhood of unity. This has a 
strong regularizing effect on the nature of the dynam- 
ics. See the surface of section in Fig. @ for the case 
(A = —0.8, e = 1). Only its spectrum gives number vari- 
ance statistics close to the GUE results; see the lower 
panel of Fig. 

To summarize, we have given an example of a simple, 
continuous, dynamical system that comes close to gener- 
ating GUE statistics. It is surprisingly difficult to find 
an essentially, fully chaotic system that does so. The pit- 
fall of false time reversal breaking can lead to symmetries 
that are quite well hidden, and the addition of a vector 
potential to a dynamical system has a strong tendency 
to move the system away from fully developed chaos. 
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